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1! ; Detailed Current Edit History 
Exere / G COMPLEX*16 giving G COMPLEX*16 result 


Page 
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Vite OSSRIvCG - G Cony Grete / G COMPLEX*16 DIVISION ROUTINE 
IDENT /1-001/ ; File: OTSDIVCG.MAR 


COPYRIGHT («) 1978, 1980, 1982, 1984 BY 
DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 
ALL RIGHTS” RESERVED. 


& : 
$3 
3 @ 
§ ;* ca 
;* ® 
ge e 
h 3* « 
11 ;* THIS SOFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED * 
\¢ z* ONLY IN ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE * 
13 ;* INCLUSION OF THE ABOVE COPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTHER * 
14 ;* COPIES THEREOF may NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO AN * 
15 ;* OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY e 
00 16 i TRANSFERRED. * 
3 e 
° 18 3* THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE * 
19 ;* AND SHOULD NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT * 
0 9 :* CORPORATION. ® 
3 & 
° ; 3* DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS * 
oS ? :* SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. . 
000 5 | sadneusanbooenenenneneenbeneseenncnemsennnnneatoannensenedannnnennsonandenes 
0000 8: 
9000 9; FACILITY: MATH LIBRARY 
000 0 ;+ 
900 1; ABSTRACT: 
000 : : Perform G COMPLEX*16 division 
000 4; 
bo 5 i: 
0 6; 
000 7 ; AUTHOR: 
4 3 3 Steven B. Lionel, 12-July-1979 
$609 0 : MODIFIED BY: 
4 1 : 
$8 : : 
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$3 -SBTTL HISTORY ; Detailed Current Edit History 
is Edit History 


39 ; 1-001 = Adapted from OTS$DIVC version 1-003. SBL 12-July-1979 
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3 -SBTTL DECLARATIONS 


: ; INCLUDE FILES: 


: EXTERNAL SYMBOLS: 


. 
. 


60 
] 
3 
8 rt: ; MACROS: 
0 ee: 
oh 08 ; PSECT DECLARATIONS: 
0000000 i -PSECT _OTSSCODE PIC, USR, CON, REL, LCL, SHR, 
00 71 EXE. RD, NOWRT, LONG 
44 ts 
$09 te : EQUATED SYMBOLS 
000004 $09 fe ; ‘ l t of dividend 
3 = ifr 
000000¢ 0000 7 b = 12 S imag port of dividend 
00000014 0000 78 C = 2 : real part of divisor 
0000001C 44 79 d = 28 ; imag part of divisor 
00 0 
000 i 
000 ; 3; OWN STORAGE: 
000 3 
000 4; none 
0000 5 
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; FUNCTIONAL DESCRIPTION: 
OTSSDIVCG_R3S = G COMPLEX*16 / G COMPLEX*16 giving G COMPLEX*16 result 


The COMPLEX*16 result is computed as follows: 

1) Let (a, b) represent the COMPLEX*16 dividend. 

§} Let (c, d) represent the COMPLEX*16 divisor. 
>) Let (r, i) represent the COMPLEX*16 quotient. 

Then: 


r = (ac + bd) / ( 
i = (be = ad) / ( 


CALLING SEQUENCE: 

Complex_quotient.wgc.w = OTSSDIVCG_R3(dividend.rgc.v, divisor.rgc.v) 
INPUT PARAMETERS: 

Dividend and divisor parameters are represented as 

FORTRAN G COMPLEX*16 numbers anG are CALL BY VALUE. 

Passing 128 bit quantities by value is a violation 

of the VAX calling standard, but is excused because 

this is a code support routine not meant to be 

callable by users. 


IMPLICIT INPUTS: 
NONE 
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OUTPUT PARAMETERS: 
NONE 

IMPLICIT OUTPUTS: 
NONE 


FUNCTIONAL VALUE: 


The G COMPLEX*16 yplus returneG is (a, b) / (c, d) 
in resisters RO-R3! This is a violation of the VAX 
calling standard, but is excused because this is 

: code support routine, not meant to be callable 

y users. 


SIDE EFFECTS: 
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Modifies registers RO-R3! , 

SS$_ROPRAND if either argument is a reserved operand. 
SSS_FLTOVF if froeting over flow 

SS$S_FLTDIV if divide by zero 
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bc 11:37:37 Umtanne.saedorspivec.mar;1 “28° (3) 


3, “M<R4,R5,R6,R7,R8,R9I,R10,R11> 


16 giving G COM 


ENTRY OTSSDIVCG 
MTHSFLAG_JACKET 


G*MTHSSJACKET_HND, (FP) 


; Perform scaling of all operands before division 


51 
50 


2$: 


owor@qa-aomenn 
LFomnmoLr nn 


WAWNTOIPINININ | | “DOCOOOOCOOOCOCOOo 


VIMO DLOOCDBLOMWWIMMWOVIOPAS LOS POuInmovoovowvnoo 


SDooooooooooooo 


VOAAAANMNMNNIT  -VIVIVIVI MOVIMIOUOO 
OOPrPOLOranrwrf GOOCO —Lf£0—WOLS 


PAPA MAVIVIUIVIUIV 
MOOY ONRWLOPOO 
Fa at at at ot ot at ee ot 
SROLLSVLLOWWw www 
TTT HHH i Dee Be Be | 
cvVvVvV"0000"00 VUUO 


o 
Pd 
a at st 2 ss — 2 2d sd 2 — s ss —) —) 2 2 ss 


SSIS NSS SA AAA AAA AAA MMIII 
NAUSEA 9 OO NOU EWN 0 OO NOAUE WN OOONO 


SooQoooooooooooo 


; establish math error handler 


; set handler address to jacket 
3; handler 


R1 = c(AP)<exp> 0, 
RO = d(AP)<exp> 0, 
3 RO = MAX (c<exp>, d 


Teseokt it 

| Ry 

<exp>) 

RO = scaling exponent 0,3777,3776,...,1 
build a floating scale factor 


scale all operands 
RB-R9 gets C 
gets d 


3: R10-R1 

; a 

3 } 

; RO = ac 

3; R7 = bd 

3; RO = ac+bd 

: R4 = ad 

3; R2 = 

3; R2 = be - ad 

; RB = 

; R10 = dd 

3; RB = cc + dd 

3; RO = (ac+bd) / (cc+dd) 
3 R2 = (bew-ad) / (cc+dd) 


(RO-R1, R2-R3) = (r, i) 


8 7 
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see 138% 91:33:43 AX/VMS Macro V04-00 Page g 
Symbol table 66-SEP-1984 11:27:3 MTHRTL.SRCJOTSDIVCG.MAR; 1 (5) 
A = 
: = Soo0000¢ 
C = 14 
i") = 1c 
MTIHSSJACKET_HND eeeeeere x 3 
OTSSDIVCG_R3 00000000 RG 1 
¢ meme ans aa aaeeee £ 
! Psect synopsis ! 
¢oeesennr omen secce oe 
PSECT name Allocation PSECT No. Attributes 
- ABS. 00000000 ( 0.) 00 ¢ 0.) NOPIC USR CON ABS LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
-OTSSCODE 00000060 (¢ 109.) O01 ¢ 1.) PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 


te eocwrecas er seon eee ee as at am 


! Performance indicators ! 


deme ence een ene ser ew cecanes b 


Phase Page faults CPU Time Elapsed Time 
Initialization 33 00:00:00 .08 90:00: 00.69 
Command processing 105 9s 4 00:00:03.82 
Pass 78 00: 8708-28 00:00:02.95 
Symbol table sort ! 00:00: 8-99 00:00:00.00 
Pass 2 4 ee fe Se 00:00:01.85 
Symbol table output 00:00:00.01 00:00:00.01 
Psect synopsis output Bs So 4 SS ae Bae S 
Cross-reference output 00:00:00.00 00:00:00.0 
Assembler run totals 26 00:00:02.12 00:00:09.41 


The working set Limit was 750 pages. 

3071 bytes (6 pages) of virtual memory were used to buffer the intermediate code. 

There were 10 pages of symbol table space allocated to hold 6 non-local and 1 local symbols. 
237 source Lines were read in Pass 1, producing 11 object records in Pass 2. 

1 page of virtual memory was used to define 1 macro. 


$ewewwwwwwwn ooo oe ooo eneooonn+ 
! Macro Library statistics ; 
ewe eeoeeoeceoeseecoceescoces 


Macro library name Macros defined 


O GETS were required to define 0 macros. 
There were no errors, warnings or information messages. 
MACRO/ENABLE=SUPPRESSION/DISABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:OTSDIVCG/OBJ=OBJ$:OTSDIVCG MSRC$:MTHJACKET/UPDATE=(ENHS :MTHJACKET) +MSRC 
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